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I.  INTRODUCTION 


This  paper  characterizes  the  solutions  of  the  equation 


(1) 


which  arises  in  several  areas  of  investigation,  e.  g.,  determining  the 
maximum  range  of  an  artillary  trajectory  for  fixed  time  under  a given 
aerodynamic  drag  [1];  determining  the  residual  velocity  of  a bullet 
traveling  through  a liquid  as  a function  of  fuze  delay  time  and  striking 
velocity  [2] ; calculation  of  the  burning  rate  for  constant  frequency  in 
a T-Bumer  experiment  [3] ; and  in  the  problem  of  intense  surface  heating 
of  a slab  [4] . 


II.  CHARACTERIZATION  OF  SOLUTIONS 


We  will  now  state  and  prove  a series  of  lemmas,  which  will 

Ajx 

characterize  the  solutions  of  e = x as  a function  of  A,  The  first 
three  lemmas  will  give  the  algorithms  for  generating  the  solutions.  In 

Ax 

the  sequel  a value  of  x which  satisfies  e = x will  be  referred  to  as  a 
solution.  Also,  log(x)  will  denote  the  natural  logarithm  of  x. 

- 1 Ay . 

Lemma  1.  For  0 < A < e , the  sequence  y^+1  = e /i,  yQ  = e,  is  a 
bounded  monotone  decreasing  sequence  and  its  limit  is  a solution. 

Proof:  y^  = e > e^O  = y^.  Assume  y^  > y^  > . . . > y^.  Then 


McCoy,  H.,  Private  consultation. 

2 

Walther,  R. , et  al.  Forthcoming  report  concerning  vulnerability 
analysis  of  MiG-21  aircraft  to  20  mm  and  improved  20  mm  projectiles 
with  various  types  fuzing. 

3 

Ibiricu,  M.  M. , Forthcoming  report  on  analysis  of  T-Bumer  results 
based  on  A-13  propellant  data. 

^ Masters,  J.  I.,  "Problems  of  Intense  Surface  Heating  of  a Slab 
Accompanied  by  Change  Phase,"  Journal  of  Applied  Physics,  Vol.  21, 
1956,  pp.  477-484. 
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y^+1  = e^i  < e^i-1  = y^.  Hence,  by  induction  the  sequence  is  monotone 

' • ’ ' • A 

decreasing  and  bounded  by  zero;  in  fact,  it  is  bounded  below  by  e . 

Therefore,  y^  +y  and  y is  a solution. 

Lemma  2.  For  0 < A < e \ the  sequence  y^+1  = A~*log(y^) , yQ  = A 
is  a bounded  monotone  increasing  sequence  and  its  limit  is  a solution. 

Proof:  For  0 < A < e"1,  y^  = A_1logCA~^)  > A~^log(e)  = A \ = yQ.  . 

Assume  yQ  < y1  < ...  < y±,  then  y.+1  = A“1log(yi)  > A^logCy^,  1)  .=  y 
Thus,  by  induction  the  sequence  is  monotone  increasing. 

To  show  the  sequence  is  bounded,  we  will  generate  a bound  for  the 
first  three  terms,  assume  a form  of  the  bound  on  y.  , k = 1,  . ..,  i,  then 

-1 

use  induction.  To  simplify  notation,  let  U = log(A  ).  Then 

X1  = A~y 

y2  = A_1logCy1)  - A-1 [U  + log (U)] , 

y3  = A_1logCy2)  = A_ 1 [U  + log (U  + log(U))] 

= A-1 [U  + logCU)  + log (1  + U_1logCU))]. 

Since  log(l  + x)  < x for  x > 0,  then 

y3  < A'l[U  + log(U) [1  + if1]]. 

Note  that  U > 1,  this  implies  log(U)  > 0. 

Assume  the  form  of  the  bound  holds  for  y. , that  is, 

yi  < A1 [U  + log(U) [1  + if1  + ...  + U2~X]]. 

Then 

yi+1  = A‘1log(yi) 

< A-1 [log(A-1[U  + logCU)[l  + U'1  + ...  + U2-1]])] 
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< A-1 [U  + log(U  + log(U) [1  + U"1  + ...  + U2'1])] 

< A-1 [U  + log(U)  + log ( 1 + log(U) [U-1  + ...  + U1-1])] 

< A_1  [U  + log(U)  [1  + if1  + ...  + U1"1]]. 


Thus,  the  same  form  of  the  bound  holds  for  y . .. 

1+  1 


Now 


E”"k  < Eu 


-k 


k=0 


k=0 


U - 1 * 


then  for  all  i 


y.  < A‘J[U  + U(U  - 1)  "*1log(U)  ] . 

> 1 Av 

Therefore,  -*■  y which  satisfies  A log(y)  or  y = e 7 , that  is,  y is  a 
solution. 

Ay 

Lemma  3.  For  A < 0,  the  sequence  y^+j  = e 7i,  y^  = e,  has  a unique 
limit  and  its  limit  is  a solution. 


Proof:  By  induction  we  will  show  the  terms  of  even  subscripts  form 

a bounded  monotone  decreasing  subsequence  and  the  terms  of  odd  subscripts 
form  a bounded  monotone  increasing  subsequence.  Now 

A p 

yx  = e < e = yQ 
y2  - eAyl  > eAy0  = yl , 


that  is, 


yl  < y2  < 6 = V 

and 

Xl  = e%  < ^2  - y3  = f&Z  < eA^l  - yr 
Assume  for  i even 


r.  . < y.  < y, 
l-l 


i-2  * 
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Then 


y._l=  exp(Ayi2)  < expCAy^  = yi+1  = exp(Ayi)  < exp (Ay .^j)  = y±. 
For  i odd  assume 


Then 

yi  = exp(Ayi  l)  < expCAy^  = y.  + 1 - expCAy.^  < exp(Ayi  2)  = yi_J. 
Hence,  by  induction  the  sequence  has  the  following  ordering: 


yl  " X3  " y5  " " y2i+l 

Suppose  y2i+1  -*■  yQ  and  y n + yE> 
Rolle's  theorem,  there  exists  x^ 

This  implies  exp(Ax^)  = A ^ < 0. 


< ...  < y2i  < ...  < y4  < y2  < y0. 

where  y^  and  y^  are  distinct.  Then  by 
£ (y0,yE)  such  that  - x)lx=Xi  , „ 

This  is  a contradiction. 


Lemma  4.  For  A < 0 there  exists  a unique  solution. 

Proof:  This  follows  immediately  from  lemma  3 and  Rolle's  theorem. 

Lemma  5.  For  0 < A < e * there  exist  exactly  two  solutions . 

Proof:  From  lemma  1 and  lemma  2 there  exist  two  distinct 

Ax 

solutions,  call  them  x^  and  x2.  Let  h(x)  = e - x.  Suppose  there 

dh  Ax 

exists  a solution  x_  distinct  from  x,  and  x„.  Now  -r—  = Ae  - 1,  that 

3 1 2 dx 

dh  _ i _ ^ 

is,  ^-vanishes  only  at  x = A_1log(A  ).  This  is  a contradiction,  since 

dh 

Rolle's  theorem  says  must  vanish  for  at  least  two  distinct  points. 


Lemma  6.  For  A > e 1 there  exist  no  solutions. 

Ax 

Proof:  Let  g(x)  = e and  f(x)  = x.  Now  g(x)  > 1 + g'(0)x  and 
g'(0)  > e *.  Then  for  a point,  xQ>  of  intersection  (a  solution)  of  f(x) 
and  g(x),  xQ  > 1.58.  For  x > 1.58,  g(x)  > g(1.58)  + g'(1.58)(x  - 1.58). 
This  implies  xQ  > 4.86  > e.  Since  g(e)  > e and  since,  for  x > e. 
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g' (x)  > g'(e)  > 1,  there  exist  no  solutions. 

Lemma  7.  For  A = 0 or  A = e , e =x  has  a unique  solution. 

Proof:  Trivially,  A = 0 has  the  unique  solution  x = 1.  For 
A = e \ x = e is  a solution.  By  assuming  another  solution  x^  ? e,  this 
results  in  an  immediate  contradiction  to  Rolle's  theorem. 


III.  CONCLUSION 


The  previous  seven  lemmas  have  characterized  the  solutions  of 
Ax 

e = x.  Figure  1 gives  a graph  of  the  solutions,  and  a tabular  form 
of  the  solutions  is  given  in  Table  1. 


There  are  more  general  forms  of  equations  for  which  solutions  are 
desired.  By  an  appropriate  change  of  variable,  these  can  be  readily 

Ax 

reduced  to  the  form  e = x.  For  example,  consider  the  solutions  of 


x + 


a = 


be 


c(x+d) 


(2) 


Let  z = ^ ^0ec^a~d^  or  x = ec^d-a^bz  - a,  then  (2)  reduces  to 

z = eAz,  where  A = bce0^-3^.  Another  form  is 


ax  = exp (bx  ) , 


(3) 


then  x = a *z^C  ^ , where  z is  a solution  of  eAz  = z for  A = bca  C, 


For  A < 0,  lemma  4 shows  the  existence  of  a unique  solution. 
Lemma  3 gives  a simple  algorithm  for  approximating  the  solution,  but 
for  A < -1,  this  algorithm  converges  slowly.  Fritsch,  Shafer  and 
Crowley  [5]  considered  the  case  for  A < 0 and  gave  an  efficient 
algorithm  for  approximating  the  solution. 


Fritsch,  F.  N.,  Shafer,  R.  E.  and  Crowly,  W.  P.,  "Solution  of 

Transcendental  Equation  weW  = x".  Communication  of  the  ACM,  Vol . 16, 
1973,  pp.  123-124. 
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TABLE  1.  SOLUTIONS  OF  THE  EQUATION  X = e 


AX 


0 . 365 

2.41047 

3. 097  68 

0.360 

2. 23913 

3. 39658 

0.355 

2.  13025 

3. 637  30 

0 . 350 

2.  04754 

3.85633 

0.345 

1.98002 

4.0  6486 

0.340 

1.92263 

4. 26824 

0.335 

1 .87257 

4.46944 

0.330 

1.82808 

4.67  0 51 

0 . 325 

1 .7880  1 

4.87  287 

0 . 320 

1 .75153 

5.07765 

0.3IS 

1.7  1804 

5 . 28  57  7 

0.310 

1 . 68707 

5. 49803 

0. 305 

1. 65826 

5.7  1515 

0.300 

1.63134 

5.93779 

0.295 

1 . 60607 

6.16660 

0 . 29  0 

1 . 58226 

6.40223 

0.285 

1 . 5597  6 

6. 64529 

0 . 280 

1 . 53842 

6*89645 

0 . 27  5 

1.51815 

7. 15638 

0.270 

1.49883 

7 . 4 257  6 

0.265 

1.48039 

7.70533 

0.260 

1 .46274 

7.99586 

0 . 255 

1 . 44584 

8. 298  18 

0. 250 

1 . 42961 

8.61317 

0 . 245 

1 .41402 

8.94177 

0.240 

1 . 3990 1 

9 . 28500 

0 . 235 

1 . 38454 

9 . 64397 

0 . 230 

1 . 37  0 58 

10.0 1987 

0 .255 

1 . 357  1 0 

10.41401 

0 . 220 

1 . 34406 

1 0 . 8 £7  6 1 

0.215 

1 .33144 

1 1 . £6283 

0.210 

1.31921 

1 1 .7  2078 

0 . 205 

1 . 307  36 

12.20354 

0. 200 

1 . 29586 

12.7 1321 

0.195 

1 . 28469 

1 3. 25207 

0 . 190 

1 . 27  38  3 

13.82269 

0.185 

1 .2  6327 

1 4 . 4 27  9 2 

0.180 

1 . 25300 

1 5 . C 7 0 9 3 

C . 175 

1 .24  30  0 

15.7  5529 

0.170 

1.23325 

1 6.48  5 0 1 

0 . 1 65 

1 . 22  37  5 

17.26459 
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TABLE  1.  SOLUTIONS  OF  THE  EQUATION  X = e (CONTINUED) 


A 

X, 

x„ 

1 

2 

0.160 

1.21448 

18-09916 

0.155 

1.20544 

18.99452 

0.  150 

1 . 19661 

19.9  57  30 

0.  145 

1. 18798 

20.99510 

0 . 140 

1 . 17956 

22.11665 

0.135 

1 . 17  132 

23.33206 

0.130 

1. 16326 

24.6531 0 

0.125 

1.  15537 

26.09349 

0.120 

1.147  65 

27. 669  39 

0.115 

1 . 14009 

29. 3999  3 

0.110 

1 . 1 3269 

31. 30792 

0.105 

1 . 1 2544 

33.42073 

0.100 

1 . 1 1833 

• 35.77 1 52 

0.099 

1.11 692  . 

36.  27  365 

0 • 098 

1 . I 1552 

36.787  £6 

0 . 097 

1.11413 

37 . 3 1 £7  3 

0.096 

1 . 1 1 £7  4 

37 .8  50  44 

0.095 

1.11135 

38.40083 

0.094 

1 . 1 0998 

38.96432 

0.093 

1 . 1 0860 

39 . 54 1 37 

0.092 

1 . 1 0723 

40 . 1 3245 

0.091 

1 . 10587 

40.73805 

0.090 

1 . 10452 

41 . 35870 

0.089 

1 . 10316 

41 .99493 

0.088 

1 . 10  182 

42.  647  32 

0 . 087 

1 . 1 0047 

4 3.  31  647 

0*086 

1.09914 

44.00300 

0.085 

1.09?  6 1 

44.7  07  56 

0 . 084 

1 .09646 

45. 43085 

0.063 

1.09516- 

46.  17359 

0.088 

1 . 09  384 

46.93654 

0.081 

1.09253 

47 .72051 

0.080 

1 . 09  I 2 2 

48. 52633 

0 . 079 

1 . 06992 

49.35490 

0 . 078 

1.08862 

50. 2 07  1 5 

0.077 

1.08733  . 

51  - 0840  6 

0 . 076 

1 .08604 

51.98668 

0.075 

1 . 0847  6 

52. 9 1 61  0 

0.074 

1 . 08348 

53. 87349 

0.073 

1 . 08221 

54.860  08 

0.072 

1 . 08094 

55.877  16 
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TABLE  1.  SOLUTIONS  OF  THE  EQUATION  X = e (CONTINUED) 


A 

xi 

X2 

0.071 

1 .07967 

56.92612 

0.070 

1 . 07641 

58.00840 

0.0  69 

1.077  16' 

.59.1 2555 

Q.Q68 

1 . 07590 

60 . 27923 

D.067 

1 . 07466 

61 .47  1 17 

0.066 

1 . 07342 

62 .7  0323 

0.065 

1 . 07216 

63.97738 

0 .0  64 

1 . 07094 

65.  29  57  3 

0.063 

1 . 0 697  1 

66. 66052 

0 . 06? 

1 . 0 6349 

68. 0741 6 

0.061 

1 . 0 67  27 

69 . 5 39  1 9 

0.0  60 

1 . 00605 

7 1.058  o'C 

0.059 

1 . 0 6464 

7 2.  634t  0 

0 . G 56 

1.06363 

74. 27  1 06 

0.057 

1 • 0 624  3 

7 5.97 1 1 1 

0 . osc 

1.06123 

77  • 7 38  38 

0.0  55 

1.06004 

7 9. 57  677 

0 .0  54 

' 1 . 0 5834 

31.49049 

0.053 

1 . 0 5766 

83.48407 

0 . 0 52 

1 . 0 5647 

35*  5624  3 

0.051 

1 .05529 

87-7  308^ 

0 .050 

1.05412 

89.9951 1 

0.049 

1 .05295 

92.36141 

0.045 

1.05175 

94.83656 

0.0  47 

1.05062 

97 .42794 

0.046 

1 . 04946 

100. 14359 

0.045 

1 . 04830 

102.99232 

0.044 

1 . 04715 

1 05.98377 

0 . 043 

1.04601  ■ 

109.1 2852 

0.042 

1 . 04486 

1 12.43818 

0.041 

1.04372 

1 15.92557 

0.040 

l .04259 

1 19.60483 

0.0  39 

1.04145 

123*49162 

0.036 

1 .040  32 

1 27. 60  333 

0.0  37 

1 .039  20 

131.959  28 

0.0  36 

1.03808 

136.58106 

0.035 

1 .03696 

141.49283 

0.0  34 

1 . 0 35S  5 

1 46. 7 21  69 

0.033 

1 . 0347  4 

1 52. 298 19 

0.032 

1 .03363 

158. 25686 

0.031 

1 . 03253 

1 64.63684 
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TABLE  1.  SOLUTIONS  OF  THE  EQUATION  X = eAX  (CONTINUED) 


A 

X1 

X2 

0.030 

1.03143 

17  1.48276 

0 . 0 29 

1.030  33 

178.84562 

0 • 0 28 

1 . 029  24 

186.78404 

0.  0 27 

1 . 0 28  1 5 

195. 36567 

0.026 

1 . 0 27  0 6 

204.66900 

0.0  25 

1 .02598 

214.78562 

0.0  24 

1 .02490 

225.82297 

0.023 

1 .02383 

237.90800 

0.0  22 

1 .02276 

251.19164 

0.021 

1 . 0 2 1 69 

265.85477 

0.0  20 

1.02062 

282.11590 

0.019 

1 .0  1956 

300.24141 

0.0  16 

1.  0 1850 

320.559  29 

0.0  17 

1 . 0 1745 

343.47779 

0.016 

1.01640 

369.51134 

0.0  15 

1 . 0 1535 

399.31705 

0.014 

1 . 01430 

433.747  29 

0.013 

1.01326 

473.92721 

0.0  12 

1 . 0 1222 

521.37 198 

0.0  11 

1.011  19 

578. 1697  8 

0.0  10 

1.01015 

647 .277  51 

0.0  09 

1.00912 

733.0  1908 

0.008 

1.00810 

641 .9677  1 

0.007 

1 .00707 

984.60593 

0.006 

1.00605 

1 178.69368 

0.005 

1.00504 

1456.79943 

0.004 

1 .00402 

1885.4851  1 

0.003 

1.0030 1 

2624.  17  38  1 

0.002 

1 .0020  1 

41  67. 5407  3 

0.001 

1.00100 

91  18.00663 

0.00 

1.00000 

-.01 

0.99015 

-.02 

0 .98058 

-.03 

0.97 128 

- . 04 

0 .96224 

-.05 

0.95345 

-.06 

0.94488 

-.07 

0.93654 

-.08 

0 .92842 

-.09 

0.9  20  49 

-..10 

. 0.9127? 

•- 
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TABLE  1.  SOLUTIONS  OF  THE  EQUATION  X = 


A X: 


1 1 

0.90522 

12 

0.89786 

-.  13 

0.89067 

14 

0.88364 

-.15  , 

0. 87  677 

16 

0.87005 

-.17 

0 .86347 

-.18 

0.85704 

-.  19 

0 .85075 

-.20 

0.84458 

-.21 

0 .83854 

-.22 

0.83262 

-.23 

0.82682 

24 

0.821 13 

-.25 

0.81555 

-.26 

0.81008 

-.27 

0.80471 

-.28 

0.79944 

-.29 

0.79427 

-.30 

0.78918 

-.31 

0.78419 

-.  32 

0.77929 

-.33 

0 .77447 

-.34 

0.76973 

-.35 

0 .7  6508 

-.36 

0 .76050 

-.37 

0 .75600 

-.38 

0 . 7 5 1 57 

-.  39 

0.74721 

- . 40 

0 .74292 

-.41 

0.73870 

-.42 

0.73454 

-.43 

0.73045 

-.44 

0.72642 

-.45 

0.72245 

-.46 

0.71854 

-.47 

0.7 1469 

-.48 

0.71090 

-.49 

0 .707  15 

-.50 

0.70347 

-.51 

0.69983 

X 
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52 

0 . 69625 

53 

0.69  27  1 

54 

0.689  23 

55 

0 . 68  57  9 

56 

0. 68240 

57 

0.67905 

58 

0. 67  57  5 

59 

0. 67  249 

60 

0. 66927 

61 

0 . 66610 

62 

0. 66296 

63 

0 . 65987 

64 

0 . 65681 

65 

0 . C5379 

66 

0 • 65081 

67 

0 . 647  87 

68 

0. 64496 

69 

0 . 64208 

70 

0 . 639  24 

7 1 

0.63644 

72 

0.63366 

73 

0. 63092 

74 

0 . 62821 

75 

0 . 62553 

76 

0.62289 

77 

0 . 62027 

78 

0 .61768 

79 

0.61512 

80 

0 . 61259 

81 

0.61008 

82 

0 . 607  60 

83 

0.60515 

84 

0.60  27  3 

85 

0.60033 

86 

0 • 59795 

87 

0 . 59561 

88 

0 . 59328 

89 

0.59093 

90 

0.58870 

'91 

0 . 58  645 

92 

0.58422 

16 


TABLE  1 


SOLUTIONS  OF  THE  EQUATION  X « eAX 


-.93 

0. 5820 1 

-.94 

0 . 57982 

-.95 

0 . 57  7 66 

-.96 

0 . 57551 

-.97 

0 .57  339 

-.98 

0 . 57  1 29 

-.99 

0 . 5 69  2 1 

-1.00 

0 . 5 67  1 4 

-1.50 

0.48391 

-2.00 

0.42630 

-2.50 

0.38343 

-3.00 

0.34997 

-3.50 

0.32294 

-4.  00 

0. 30054 

-4.50 

0 . 28  1 61 

-5.00 

0 .26534 

-5.50 

0.251 19 

-6.00 

0. 2387  3 

-6.50 

0.227  67 

-7.00 

0. 21776 

-7.50 

0.20883 

-8.00 

0.20073 

-8.50 

0.  19333 

-9.00 

0 • 18656 

-9.50 

0 . 18032 

-10.00 

0.17455 

-10.50 

0 . 169  20 

-11.00 

0.1 6423 

-11.50 

0.15958 

-12.00 

0. 15523 

-1 2.50 

0.15116 

-13.00 

0.  147  32 

-13.50 

0.  1437  0 

-14. 00 

0.  140  29 

-14. 50 

0. 13706 

-1 5. 00 

0 . 1 3400 

-15.50 

0 . 13109 

-16.00 

0 .1 2832 

-16.50 

0. 12569 

-17. -0  0 

0.12318 

-17.50 

0 . 12079 
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TABLE  1.  SOLUTIONS  OF  THE  EQUATION  X = e (CONTINUED) 


A 

X1 

-18.00 

0. 1 1849 

-18.50 

0 . 1 1630 

-19.00 

0. 1 1420 

-19.50 

D . 1 1218 

-20.00 

0. 1 1025 

-20 . 50 

0.  10839 

-21.00 

0 . 10660 

-21.50 

0. 10488 

-22.00 

0 .10322 

-22. 50 

0.10162 

-23.00 

0.10008 

-23.50 

0 . 09859 

-24. 0 0 

0.09715 

-24. 50 

0 . 09575 

-25.00 

0.09441 

-25.  50 

0.09310 

-26. 00 

0.09184 

-26.  50 

0.09061 

-27.00 

0 . 08942 

-27. 50 

0*08827 

-28.00 

0.0  87  1 5 

-28.50 

0 . 08606 

-29.00 

0.08500 

-29. 50 

0.08397 

-30.00 

0.  0 8 297 

-30.50 

0 . 08200 

-31.00 

0.08105 

-31.50 

0.08013 

-32.00 

0.07923 

-32.  50 

0.07835 

-33.00 

0.  07750 

-33.  50 

0. 0 7 667 

-34.00 

0 .07  585 

-34.50 

0 .07506 

-35.00 

0.07428 

-35.  50 

0.07352 
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